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ABSTRACT 

Wc develop one-zone galaxy formation models in the early Universe, taking into account dust 
formation and evolution by supernova (SN) explosions. Especially we focus on the time evolution 
of dust size distribution, because H2 formation on the dust surface plays a critical role in the 
star formation process in the early Universe. In the model we assume that star formation rate 
(SFR) is proportional to the total amount of H2. We consistently treat (i) the formation and 
size evolution of dust, (ii) the chemical reaction networks including H2 formation both on the 
surface of dust and in gas phase, and (iii) the SFR in the model. First, we find that, because of 
dust destruction due to both reverse and forward shocks driven by SNe, H 2 formation is more 
suppressed than that without dust destruction. At the galaxy age of ~ 0.8 Gyr, for galaxy models 
with virial mass M vir = 10 9 Mq and formation redshift z v - u — 10, the molecular fraction is 2.5 
orders of magnitude less in the model with dust destruction by both shocks than that in the 
model without dust destruction. Second, wc show that the H2 formation rate strongly depends 
on the ISM density around SN progenitors. The SFR in higher ISM density is lower, since dust 
destruction by reverse shocks is more effective in higher ISM density. We conclude that not only 
the amount but also the size distribution of dust being related with the star formation activity 
strongly affects the evolution of galaxies in the early Universe. 

Subject headings: dust — galaxies: evolution — galaxies: formation — galaxies: ISM — early Universe 
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1. Introduction 

Understanding of galaxy evolution in the early 
Universe remains one of the most important goals 
of modern astrophysics. Modeling of primeval 
galaxy formation requires an accurate treat- 
ment of star formation p rocess in low-metallicity 

Glover fe Jappsenll2007t 



gas (IJappsen et al. l200Yt juiover &l jappL- 
Smith et al.l 120081 l2009t IJappsen et al.l l2009allbl ) 



A critical challenge for achieving this goal is 
due to our poor understanding of how gas is 
converted into stars under different conditions 



dKrumholz fe McKeel2005HRobertson fe Kravtsov 
l2008t iGnedin et al.ll2009l ). In particular, star for- 
mation efficiency in primeval galaxy is still uncer- 
tain. 

The standard approach in theoretical stud- 
ies of galaxy formation so far is to adopt a 
recipe which ties the star formation rate (SFR) 
to gas den s ity bo th in semi-analytic models (e.g. 
ICole et all l2000h a nd in numerical si nulations 



(e.g. iSpringel et al.ll2005l) . Such a recipe is based 
on the empirical correlations observed in lo- 
cal galaxies, nam ely the Kcnnicut- Schmidt law 
(|Kennicuttlll998l ). These correlations have only 
been studied relatively well for nearby massive or 
star bursting galaxies. However, for galaxies with 
low surface brightness and/or low-metallicity, this 
empirical relation may not be yalid. Indeed, both 



nearby metal-poor galaxie s dBigiel et al 



2008) 



and high-redshift galaxies (jWolfe fe Chen! |2006[ ) 
provide a variety of clues suggesting that gas 
conversion into stars in low-mass, low-metallicity 
galaxies is very inefficient. 

The star formation efficiency may depend on 
ability to convert a fraction of gas mass into molec- 
ular form. Molecular hydrogen is produced by 
chemical reactions in gas phase in first galaxy 
halos. In the reionization era, H2 molecule dis- 
sociation by the Lyman- Werner ultraviolet (UV) 
background between 11.2 and 13.6 eV is impor- 
tant in the lower mass H2 cooling halos. Gas 
condensation in the lower mass H2 cooling halos 
can be delayed by the Lyman- Werner background 



(iMachacek et al.ll200ll l2003t lYoshida et alil2003t 



200 



Susall2007t IWise fc Abel 20071 lO'Shea fe Norman 



The Lyman- Werner background thus in- 
creases cooling times in the centers of such ha- 



1 yamasawa® astrol.sci.hokudai.ac.jp 



los. As a result, the minimum mass of a star- 
forming halo increases with the Lyman- Werner 
background intensity. The Lyman- Werner back- 
ground becomes less of an issue in atomic line cool- 
ing halos as Lya cooling provides ample amounts 
of free electrons for H2 coolin g, and they become 



2009). 



self- shielding to this radiati on fO'Sh ea fe Norman 
20081: ISusall2008t I Wise fc Abell 120081: 1 Wise fc Gen 



In the later epoch, dust ejected by stars in 
galaxies is effective to shield the Lyman- Werner 
background and acts as an effective catalyst for 
H2 molecule production on the dust grains. In 
simulations with sta r formation models based on 
molecular hydrogen ( Robertson fc Kravtsovll2008t 



Gnedin et al. I l2009h . once the gas enriched up to 



Z^O.01 — 0.1 Zq, the subsequent star formation 
and enrichment o f metal and dust ca n be much 
more accelerated. Gnedin et al. I (120091) show that 
the transition from atomic to molecular hydrogen 
depends primarily on metallicity, assuming that 
the dust abundance is directly related to metallic- 
ity. 

Dust plays a crucial role in the star forma- 
tion: (i) molecular hydrogen is produced more 
efficiently on dust grains than in gas phase, (ii) 
dust shields dissociating UV radiation, and (iii) 
dust allows the formation of low- mass stars in low- 
metallicity environments, an d hence affects the 



initial mass function (IMF) (Omuk ai et al 



Schneider et all l2006t ISchneider fc Omukai 
Omukai et al.ll2010h . 



2005 



2010: 



In theoretical studies on the molecular abun- 
dance in the interstellar medium (ISM), dust 
abundance is often scaled with the metallicity and 
dust grain properties are assumed to be the same 
as in the local ISM. However, the composition 
of dust is likely to be different in early galaxies. 
The observational evidence is that the dust extinc- 
tion curves of the broad absorption line quasars 
at z > 4 are likely to be due to the type II SN 
(SN I I) dust (|Maiolino et alj|2004 iGallerani et aT 
2010h . 



Since the lifetime of SN II progenitor is short, 
SN II can be the dominant production source 
of dust grains in young (< 1 Gyr) galaxies. 
Primeval SNe produced by Population III stars 
(iBromm et al l (20031 iKitavama fc Yoshidal 12005 ; 
IWhalen et all l2008h may contribute the dust 



production (|Nctzawa et ah 2003 : Schneider et al 
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2004j). The winds of evolved low-mass stars con- 
tribute to dust formation considerably in nearby 
galaxies, but the cosmic time is not long enough 
for such stars to evolve at high redshift (z > 5) 
where all galaxies should have ages younger than 
~ 1 Gyr. Contribution of dust production by low- 
mass stars is not dominant in such young galaxies. 
In addition, dust is destroyed by SN shocks. Thus, 
the modeling of dust evolution in galaxies requires 
an accurate treatment of production and destruc- 
tion of d ust grains together with st ar formation 
activities ( Hirashita fc Ferrara 20021 ). 

In this paper, we investigate not only the evo- 
lution of dust mass but also the time evolution 
of dust size distribution. The dust size distri- 
bution evolves rapidly because of the destruc- 
tion by sputtering in the high-velocity shocks 
driven by SNe. Collision of the expanding 
SN ejecta with the surrounding ISM creates 
a forward shock at t he interface betwee n the 
ejecta and the ISM ([Nozawa et al.l 120061 ). and 



a reverse shock that penet r ates into the ejecta 



N 



Bianchi fc Schneider! 120071: iNozawa et ail 12007 



Nath et al.ll2008USilvia et al.ll2010Q . Since the ero- 
sion rate by sputtering does not strongly depend 
on the grain size, small grains are predominantly 
destroyed regardless of grain species. Therefore, 
the fraction of small size grains relatively decreases 
with galaxy evolution. 

We focus on the effects of molecular hydro- 
gen abundance on the SFR in the early stage 
of galaxy evolution, taking into account molecu- 
lar formation on dust, s ince H2 formation on dust 
surface is very effe c tive ( Hirashita fe Ferrarall2002l: 
Cazaux fc Spaansl 20041 1 Irlirashita fe Ferrara 



(l2002h show that this effect causes an enhance 
ment of the SFR by an order of magnitude on 
a timescale of 3 — 5 galactic dynamical time. 
However, they assumed a single dust grain size 
(~ 0.03 /zm). We adopt more accurate analytic 
formulae for the form ation of molecular hydrogen 
on dust grains than Hirashita fc Ferrar 1 (l2002l) 



by using the r esults of du st size distribution by 
Nozawa et all (|2006l l2007h . 



This is the first study on galaxy evolution con- 
sidering dust size evolution for halo masses above 
10 8-9 in the high-rcdshift (5 < z < 10), whose in- 
teriors we expect to be roughly self-shielded from 
both ionizing and Lyman- Werner UV radiation. 
To show clearly the dependence of galaxy proper- 



ties on dust destruction, we use a simple one-zone 
galaxy model. 

The paper is organized as follows. In §2 we 
describe the dust evolution model. In §3 we ex- 
plain our one-zone galaxy model. In §4 we present 
the results. In §5 we discuss the effects of the 
dust size evolution on H2 formation process and 
conclude by summarizing our results. Throughout 
this paper we adopt the cosmologic al parameters 



from the third-year WMAP results (ppergel et al 
20071) . n A = 0.76, ft M = 0.24, n b = 0.04, and 



H Q = 73 km s _1 Mpc . 

2. Dust evolution model 

2.1. Source of dust in the early Universe 

SNe II are believed to be the dominant sources 
of dust at high redshift of z > 5 because of 
short lifetimes (< 10 7 yr) of th e ir massive pro 



2007; Gallerani et al 



genit ors (e.g. iDwek et al 

2010h . Dust formation in the ejecta of pri 



mordi al SNe II has been investigated theoreti- 
cally dTodini fc Ferrarall200ll : INozawa et al.lliooj 
Cherchneff fc Dwekl l2010h . The amount and 
the size distributions of dust grains injected 
into ISM have been investigated by consider- 
i ng the destruction in SN remnants (SNRs ' 



Nath et al.l 



( Bianchi fc Schneider! 12007 

20081: ISilvia et al 



Nozawa et al 
20101) 



2007 



Although 

how much dust re ally forms in the eje cta has been 
still under debate (K ozasa et al. 2009, for review) , 
the recent observations of Cas-A SNR revealed 
the pre sence of ~ 0.07 Mp dust condensed in th e 
ejecta ( Barlow et al. 201dt Sibthorpe et al. 2010f ). 
which is consistent with the dust mass predicted 
by the theoretical model taking into account for- 
mation and destruction process es of dust in a 
Type lib S N dNozawa et al.ll20ldl). IValiante et al ' 



(|2009l ) and IDwek fc Cherchnefll (|2010l ) have pro 



posed that the contribution from asymptotic giant 
branch (AGB) stars in the high-redshift quasar 
J1148+5251 cannot be neglected for the total 
dust budget even at z ~ 6. However, the size 
distribution of dust formed in the mass-loss wind 
of AGB stars has not been fully studied ye t 
(|Ferrarotti fc Gaill 120061 : IZhukovska et al.l l2008h . 
If type la SNe could occur in such an early epoch, 
they are unlikely to be efficient sources of dust 
(Nozawa et al. 2010 in preparation). Therefore, 
in order to follow the evolution of dust size distri- 
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bution and reveal the resulting influence on galaxy 
evolution, we consider SNc II as the source of dust 
in the early Universe. 

The basic quantity for governing the production 
and destruction history of dust by SNc II is the 
rate of SN II explosions, 7sn(£), given by 



7sn(*) 



dm*(t-r(m))^(m), (1) 



where is the SFR at time t, <f>(m) is the 

stellar IMF, r(m) is the lifetime of a star whose 



mass is m, and 



and 



are the upper 



and lower mass limits of SN II progenitors, re- 
spectively. In this pap er we adopt th e Salpeter 
IMF 0(m) oc to" 2 35 , (|Salpeter|[l955h ) with the 
stellar mass range between 0.1 Mq and 60 M Q , 
we assume mj, M = 8 M Q and mg N = 40 Mq 
(|Herger et all l2003h ■ For r(m), we adopt the 
model of zero-m etallicity stars without mass loss 
(ISchaererll200l . 

In this paper, we do not consider Population 
III stars, for simplicity. In the forthcoming pa- 
per, we will consider possible contribution of Pop- 
ulation III stars. Population III stars formed 
out of the primordial gas are considered to be 
much more massive than Populat i on I/ II stars 



(jYoshida et all 120081: iBromm et ail 120091 for re 



views), and thus the primordial IMF might be bi- 
ased toward a higher mass (> 10 Mq) than that 
in the present Universe. Furthermore, Population 
III stars as massive as 140 — 260 Mq are predicted 
to end their lives as pa ir-instability SNe (PISNe 



Herger k, Wooslevl 120021 ) and to produce a large 



amount of metals an d dust ( Nozawa et al.l 120031: 



Schneider et all 120041) . Ho w ever. IJoggerst et al„ 
( 2010al ) and IJoggerst et all (|2010bl ) address the 



growing nuclcosynthetic 'forensic' evidence that 
the majority of primordial stars may have been 
15 — 40 Mq objects. On the other hand, once 
the gas is enriched up to a critical mctallicity 
of Z ~ 10~ 6 — 10~ 5 Zq, formation of low-mass 
stars is triggered, leading to the transition of the 
star formation mode from massive population III 
stars to low-mass population I/II stars, if dust 
present (lOmukai et al.l2005l:ISclineider et alj200& 
ISchneider k OmukafcOloHOmukai et al.l2010l) . If 
there is no dust, the transition of the star for- 
mation mode is expe cted to occurs at 10 -35 Zq 
(jMackev et al.l 120031 ). In this case, the CMB 
limits the lower masses of stars to a few 10s 



of M Q ((Smith et al.l 120081: ISchneider fc Omukai 
2010). Formation history of galaxies considering 
the time-dependent IMF from the top-heavy to 
the Salpeter-like IMF and taking into account the 
production and destruction of dust by PISNe, will 
be explored in the forthcoming paper (Yamasawa 
et al. 2010 in preparation). 

2.2. Dust injected from SNe II into ISM 

Throughou t this paper we adopt the models by 
Nozawa et all (I2003ll2007l) for du st formation and 
destruction. iNozawa et al.l (|2003l ) investigated the 
dust production in the ejecta of primordial SNe II 
as well as PISNe, applying a theory of non-steady 
state nucleation and grain growth. They revealed 
the grain species formed in the ejecta and their size 
distributions for the unmixed and mixed elemental 
compositions within the He core. In what follows, 
we apply the results of calculation for the unmixed 
ejecta of SNe II with the progenitor mass m = 
13,20,25, and 30 Mq and the explosion energy 
10 51 erg, and extrapolate the data to the mass 
range from 8 to 40 Mq. 

The solid line in Figure [1] shows the IMF- 
averaged mass distribution of dust formed in the 
ejecta of SNe II, M®(a), which is weighted by the 
Salpeter IMF and is summed up over all the grain 
species as, 



M° d (a) = 

3 



Ej:? N dm Ml^mMm) 



E, £r dm 0(m) 



,(2) 



where M^j(a,m)da is the mass of the j-th dust 
species produced in a SN II with radii between a 
and a+da as a function of progenitor mass m, and 
superscript means the case with no dest ruct ion 
by a reverse shock. In Figure [U we plot aM®(a) 
in the vertical axis to make clear the mass fraction 
in each logarithmic bin. We can see that the grain 
radii range from a few A up to a few //m and that 
the size spectrum of dust in mass has a peak at 
a ~ 0.1 /xm. 

In the course of their injection into ISM, dust 
grains formed in the ejecta are destroyed due to 
sputtering in the hot gas between the reverse and 
forward shocks, which is hereafter referred to as 



the destruction by reverse shock. INozawa et al 
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Fig. 1. — The IMF-averaged mass distributions of 
dust formed in the ejecta and injected into ISM 
per SN II: The solid line denotes the mass distri- 
bution of dust befo re the destruction through a 
reverse shock, aA4 d (a) and the dashed lines and 
the dotted lines denote the mass distributions of 
dust at the injection into ISM after the destruc- 
tion through a reverse shock, aMj SN (a), for the 
number densities of gas around a SN progenitor, 
n SN = 0.03, 0.1, 0.3, 1.0, 3.0, and 10.0 cm- 3 , 
which are used to annotate the curves. The hori- 
zontal axis shows the radius of dust in units of cm. 
The vertical axis is the mass distribution of dust 
aMf^{a) in units of M Q . 



(|2007l ) investigated the survival of the newly 
formed dust in the shocked gas within the SNRs 
expanding into the uniform ISM with hydrogen 
number densities of nsN = 0.1, 1.0 and 10.0 cm -3 , 
and showed that the destruction efficiency of 
newly formed dust is not only sensitive to the 
initial size distribution but also strongly depends 
on tisn- To investigate the dependence of destruc- 
tion of dust on the ISM densities, we extend their 
models to six cases of n-sN = 0.03, 0.1, 0.3, 1.0, 
3.0, and 10.0 cm" 3 . 

Figure Q] shows the IMF -averaged mass distri- 
bution of dust, Mj SN (a), injected into ISM after 
destruction by the reverse shock for nsN = 0.03, 
0.1, 0.3, 1.0, 3.0, and 10.0 cm" 3 , which is weighted 
by the Salpeter IMF and is summed up over all the 
grain species as in Equation (|2|). The mass distri- 
bution of the j-th dust species for the case of nsN, 
M^ia), is the IMF-averaged one after the de- 
struction by the reverse shock. We can see that the 
change in the dust mass distribution through pro- 



cessing in SNRs becomes more (less) prominent for 
higher (lower) gas density; small size grains get de- 
ficient with increasing nsN- The dust grains with 
radii below 0.01 /*m are preferentially destroyed 
by sputtering for tisn > 0.03 cm -3 , while dust 
with radii larger than ~ 1 /mi are almost intact 
for 71sn < 3.0 cm -3 . As a result, the mass of dust 
injected into the ISM is fully dominated by grains 
with radii above ~ 0.1 /xm. 

The total geometrical cross-section of dust per 
metal mass is an important quantity for H2 for- 
mation on grain surface at a certain metallicity 
level. The total geometrical cross-section of dust 
injected into the ISM from a SN depends on the 
gas density around the SN progenitor, nsN- In 
many papers, the total geometrical cross-section 
of dust is scaled to metal mass under following two 
assumptions: (i) the depletion factor which is de- 
fined as the dust mass per metal mass, is identical 
to that in the MW, and (ii) the dust size distri- 
bution which determines dust area per unit dust 
volume (a 2 )/ (a 3 ), is the same as that in the those 
in the MW. However, the depletion factor and the 
size distribution of the SN II dust are quite differ- 
ent from those in the MW. 



Figure [2] shows the depletion factor, M^ N d /m m 
after the destruction by reverse shock as a func- 
tion of ISM density n SN , where M^ d is the IMF- 
averaged total dust mass ejected into ISM by a 
SN, 



m SN,d — 



daX" SN (a), 



(3) 



and m m is the IMF-averaged total metal mass 
ejected into ISM, 



Si J3 SN m m,i(™)0(m)drn 



(4) 



vnm.,i{m) is the mass of i-th element of metal 
ejected from SN with progenitor mass is m and 
is taken from lUmeda fc Nomotol (|2002l ). The de- 



pletion factor in the case of «sn = 1-0 cm is a 
factor 6 smaller than that in the MW. 

Figure [3] shows the ratio of total geometric 
cross-section of dust to the total metal mass, 
-Darca/WT^ after the destruction by the reverse 
shock as a function of ISM density n^. The total 
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Fig. 2.— The depletion factor, M^ N d /ra^, for 
various dust destruction models with the densi- 
ties of ISM around the SN II progenitor, nsN = 
0.03,0.1,0.3,1,3,10 cm- 3 . T he dot- dashed line 
represents the depletion factor, Mg N d /m^ for the 
model without reverse shock destruction. The 
ejected metal ma s s per SN II is taken from 
Umeda fe Nomotol (|2002l ). The dotted line rep- 
resents the typical ratio in MW, where we as- 
sume a dust-to-gas mass ratio to be V = 0.00934 
( Pollack et al.lll994h . and assume a ratio of metal 
mass to hydrogen mass to be 0.0204 (jOmukai 
2000h . 



geometrical cross-section of dust weighted by the 



Salpetcr IMF, £>arca, is written as 



t 4 * 



daM^(a)/a, 



(5) 



where pj is the bulk density of j-th dust species. 
The ratio, D^icn/rfh^, is smaller for larger n$N be- 
cause small dust grains are efficiently destroyed 
by sputtering in the SNR under large nsNj note 
that the surface area per dust mass is larger for 
smaller-size grains. Also with increasing nsN, the 
reverse shock becomes stronger and destroys dust 



by sp uttering more effectively (see iNozawa et al 
20071 ). In addition we plot the ratio for the case 



without reverse shock together with the typical 
value in the MW for Comparison. D area / TTlra in t he 
model without reverse shock is a factor of 4 smaller 
than that in the MW. D^rca/ru^ in the case with 
"sn = 1.0 enr 3 is 40 times smaller than that in 
the MW. Thus, the rescaling of cross-section of 
dust by the metal mass using the MW value re- 
sults in significant overestimate for H2 formation 



Fig. 3. — The ratio of total dust cross-section 
to total metal mass, D2rca/rrh^, for various 
dust destruction models with the densities of 
ISM around the SN II progenitor, «sn = 
0.03,0.1,0.3,1,3,10 cm' 3 . The dot-dashed line 
represents the ratio, D° rea /m m for the model with- 
out reverse shock destruction. The ejected metal 



mass per SN II is taken from Umeda fc Nomotol 
(|2002h . The dotted line represents the ratio in our 
galaxy where we assume a dust properties as a 



-3.5 



cm < a < 



size distribution , f(a) ex a 

2.5 x 10~ 5 cm) (jMathis et al.lll977l ; iDraine fc Lee 



(5x 10 



19841) and a dust-to -gas mass ratio, V = 0.00934 
( Pol lack et al I ll994l) . and assume that a ratio of 
metal mass to hydrogen mass is 0.0204 and th e 
bulk density of dust is 3.0 g cm -3 (jOmukai l2000h . 



in high-redshift galaxies. 

2.3. Destruction of interstellar dust by SN 
forward shocks 

Dust grains injected into the ISM are sub- 
jected to destruction by the blast waves (the high- 
velocity inters tellar shocks) driv en by the ambi- 



ent S Ne (e.g. I Jones et al.l 119941 ) . INozawa et al 



( 20061 ) investigated the processing of interstellar 
dust by sputtering in the hot gas swept up by the 
SN forward shock. Adopting the dust model by 



Nozawa et alj (120031) as the size distr ibution of in- 
terstellar dust, INozawa et al.l (|2006l ) have shown 



that the destruction efficiency of dust depends on 
the ISM density and the explosion energy of SNe 
as well as the initial size distribution of dust. It 
should be noted that the size distribution as well 
as the destruction efficiency changes as a function 
of time because interstellar dust are supplied from 



G 



SNe and processed in ISM successively according 
to star formation activity. Thus, we must deal 
with the destruction process in a way that is ap- 
plicable to any dust size distribution to explore 
the global evolution of dust size distribution. 

In order to evaluate the destruction efficiency 
of interstellar dust for any initial size distribu- 
tion, here we introduce the conversion efficiency 
as defined below. Consider that the j-th dust 
species residing in the ISM, whose size distribu- 
tion is given by the number of dust grains with 
radii between a and a + da, fj(a)da, is processed 
by sputtering in hot plasma produced through the 
a single passage of SN shock. The conversion ef- 
ficiency r]j{a, a') is defined as the number fraction 
of dust grains with radii between a' and a' + da' 
that are converted to grains with radii between a 
and a + da by sputtering through the passage of 
a SN shock. The number of dust grains with radii 
between a and a+da produced by the sputtering is 
given as r)j(a,a')fj(a')da'. Note that r)j(a,a') = 
for a > a' . Then the change in the number of dust 
grains with radii between a and a + da caused by 
a shock processing is given by 

oo 

dNjia) = ^tyfooO/itoOdo' 

a/>a 

- [1 - r)j(a,a)] f 3 (a)da 

i 

Vj( a , a')fj(a')da' - fj(a)da, 

(6) 

and as well the corresponding change of the mass 
is given by 



dM dj (a) 



4tt 



pj / vA a i a ')fj( a ') da ' 



(7) 



3 jo 
—Md,j{a)da 

where Mdj{a)da is the mass of the pre-shocked 
dust. The size distribution function after the 
shock processing /j (a) is given by / ■ (a) = fj (a) + 
dJVj/da. 

The conversion efficiency 77(0, a') and the mass 
of ISM gas swept up by shock M swept depend on 
the progenitor mass, expanding energy and type of 
SN as well as the structure, number density and 
metallicity of ISM gas. For these parameters of 
SNe and ambient ISM, once M swept and 77(0,0') 
for each dust species are calculated, the time evo- 
lution of dust mass and size distribution can be 



followed in a consistent way with the star forma- 
tion activity in galaxies as described in Section [3] 

The calculations of 77(0, a') and Af swcpt are 
performed by using the method developed by 
Nozawa et al.l (|2006l ) as follows; the efficiency of 
dust destruction increases with increasing the ex- 
plosion energy and / or increasing tt-sn but is almost 
independent of the progenito r mass as long as th e 
explosion energy is the same ( Nozawa et al.ll2006l) . 
We assume that SNe driving high-velocity shock 
in ISM are represented by Type II SN with the 
progenitor mass of 20 Mq and the explosion en- 
ergy of 10 51 erg. The ISM surrounding the SN 
is considered to be uniform with hydrogen num- 
ber densities tt-sn = 0.03, 0.1, 0.3, 1.0, 3.0 and 
10 cm -3 . By distributing dust grains with ra- 
dius a' uniformly in ISM, the conversion efficiency 
77(0, a') is evaluated for each grain species by cal- 
culating the erosion of dust by sputtering until the 
truncation time £tr which is defined as a time when 
the shock velocity is dece lerated below 100 km s _1 
(see Nozawa et al. 20061 for the details). In the 
calculations, the radii of grains in the ISM range 
from 0.00013 to 6.3 /jm for each grain species. 
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Fig. 4. — The dust mass, J2j a -Md,j(a), processed 
by a forward shock driven by a SN II explosion 
in the ambient ISM gas density of tisn = 0.03, 
0.1, 0.3, 1.0, 3.0, and 10.0 cm~ 3 , which are used 
to annotate the curves, and metallicity of Z = 
10~ 4 Zq. We assume the initial size distribution 
to be aMd,j{a) = 1 for each species j. In the cal- 
culations, we use the conversion efficiency 77^(0, a') 
evaluated for SN explosion energy, 10 51 erg, and 
progenitor mass, m = 20 Mq. 

In Figure 21 we present the changes in the dust 
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size distributions due to the shock-processing for 
different ISM-densitics, where the initial mass dis- 
tribution of dust is set to be aA4d,j(a) = 1 for 
clarity. As can be seen from the figure, small- 
size grains are destroyed significantly due to the 
erosion by sputtering, and more dust grains are 
processed for a higher ISM density. 

The mass of gas swept up by the forward shock 
until the truncation time tt r , M avropt , depends on 
not only the ISM density but also the initial metal- 
licity of the gas in the ISM. As the line cooling by 
heavy elements becomes more efficient for a higher 
gas metallicity, the forward shock is decelerated 
more quickly, resulting in a smaller M swept . By 
fitting M swept calculated for different tisn and Z, 
we derived the following approximation formula, 



M 



swcpt/M© = 1535n-^ 2U2 \{Z/Z Q ) + 0.039] 



-0.298 



(8) 

whose fitting accuracy is within 16% for 0.03 cm 3 < 
tisn < 30 cm- 3 and for 10~ 4 < Z/Z Q < 1.0. 

2.4. Formulation of dust size evolution 

In terms of the conversion efficiency describing 
the processing of dust by sputtering, here we for- 
mulate the time evolution of the mass of j-th dust 
grains with radii between a and a + Aa in our 
model galaxies, AMdj(a,t) = ^-a 3 pjfj(a 7 t)Aa, 
as 

dAM dd (a,t) 



dt 



AA/sN,dj(a)7SN(i) 



M, 



swept 



Mi SM {t) 



7Sn(*) x {AM dj (a,t) 



4tt 



Aa'-q 3 {a,a')f 3 {a' ,t)p J —a 3 } 



"*(*) 



AA/ d ,j(o,t) 
Misu(t) 



(9) 



where MisM.{t) is the total mass of gas and dust, 
and fj (a, t) is the size distribution function of dust 
species j in the ISM at a time t. We note that 
(M swep t7sN(i)/MisM(i))~ 1 is a timescale of sweep- 
ing whole ISM by SNc. The IMF-averaged mass 
of dust species j with radii between a and a + Aa 
injected from SNe II into the ISM with num- 
ber density njgM.sN is defined by AMsN,d.j(a) = 
(a) Aa. The first term on the right-hand side 
is the injection rate of dust from SNe II. The sec- 
ond term is the destruction rate of interstellar dust 



by SN blast waves, and the third term is the rate 
at which the interstellar dust is incorporated into 
stars. 

3. Galaxy evolution model 

3.1. Dark matter halo and physical state 
of gas 

We quantify the properties of dark matter ha- 
los, assuming a dynamically equilibrium state. 
The radius of dark matter halo, r V i r , is estimated 
in terms of the mass of dark halo, M V i r , and the 
redshift of virialization, z v j r , as 



4 3 



{i + s c (z viI )} p c0 n M (i + 2 Vir ) 3 = m v 



(10) 

where p c o = 3Hq /8ttG is the critical density of the 
Universe at z = 0, S c (z v i r ) is the overdensity of a 
dark matter halo vilialized at z v j r , and G is the 
gravitational constant. 

We assume dark matter halos as singular 
isothermal spheres and rotating uniform gas disks 
in their gravitational potentials. Cosmological 
N-body simulations show that structures of dark 
mat ter halos are well de s cribed by t he NF W pro- 
file dNavarro et al.lll996l ). iMo et all (Il998l ) studied 
a simple disk model in the gravitational poten- 
tial of the singular isothermal sphere and more 
realistic disk model in the gravitational poten- 
tial of the NFW halo profile. We adopt a radius 



of the disk, rdj s k — 0.18r v j r ( Ferrara et al. 2000t 



Hirashita fc Ferraral2002( ) , by considering the con- 
servation of angular momentum and assuming a 
typical value f or the spin paramete r A = 0.04 from 
the paper by iFerrara et al. (<2000h who estimate 



the radius of the disk as rd; s k = 4.5Ar v j r in a 
modified isothermal halo. 

In our one-zone model, we need a virial tem- 
perature for the initial gas temperature and a dy- 
namical timescale of gas in the disk. Gas collapsed 
at z v j r in the dark matter halo of M v ; r has a virial 
temperature, T v i r , defined as 



3fcB7Vir 



(11) 



where fee is the Boltzmann constant, \i is the mean 
molecular weight, and mjj is the mass of a hydro- 
gen atom. The initial value for the temperature of 
gas, T, is assumed to be T w - a . A circular velocity, 



S 



is defined as 



GAL 



(12) 



and we also define a rotation timescale, t c \ r , as 



tc 



27rrdisk 



(13) 



Note that the rotation timescale of the gas disk, 
i c i r , depend only on virialization redshift, z V u, as 



t cir = 9.3 X 10' 



yr 
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1 + z v 



18tt 2 



(14) 

We must estimate the number density of the hy- 
drogen gas, 7Jh, because it affects both the chemi- 
cal reaction rate and the cooling rate. The cooling 
time of halo gas is much shorter than the Hubble 
timescale for the objec ts of interest in this p aper 
(Tvir > 100007T) (e.g. iHutchings et all [20021) . It 
is widely understood that most z ~ 10 galaxies 
were not clear disk galaxies; numerical simulations 
that proceed from cosmological initial conditions 
(z ~ 100 — 200) clearly reveal that they possess 
highly irregular structures whose SF rates are not 
easily quantifiable, that filamentary accretion and 
frequent mergers are still churning the halo at this 
epoch, and that turbulent flows arise in the center 
of the halo that prevent cohe rent disks forming on 
the spatial scales of galaxies (I Johnson et al.ll2008 ; 



Greif et al]l20ld: IWise et al.ll20ldh . We 



assume 

that a significant fraction of baryons finally col- 
lapses to a disk in the dark matter halo potential 
for simplicity. In semi-analytic models, it is as- 
sum ed that the cooled halo gas settles into the disk 
l2000h . ar assump- 



Cole et al 



(e.g. 

tion in our model, but more detailed treatment for 
H2 formation, dust evolution and star formation 
in the gas disk. The radius of disk, ry ii sk, is deter- 
mined following iHirashita" fc Ferraral (|2002l ). We 
estimate the typic al scale height, H, from hy dro- 
static equilibrium ( Shakura fc Sunvaevlll98S 



H 



V2 — Tdisk 



2T 

3-Zvh 



^disk 



(15) 



for H/rdisk < 0.1, otherwise, assume H/r^isk = 
0.1, where v s = (fcsT//xmH) 5 is the isothermal 



sound velocity. Therefore, initial hydrogen density 
of disk, nH, is estimated as 



Mb 



(16) 



The initial mass of hydrogen in the galaxy, Mh, is 
written as 



Ah 



M s 



Ah 



m H 



' (m H + m Hc 2/Hc) 



AI v \ r flh 



&m ("Jh + m H e2/He) 



(17) 



where Mu e is the initial mass of helium in the 
galaxy, M gas = Mh + Mh , "t-Ho is the mass of 
a helium atom and yn c is t he helium abu ndance. 
We assume y H c = 0.0972 (jOmukaillioOoh . Note 
that since massive star ionizes surro unding gas and 
forms an expanding H II reg ion (e.g. Whalen et al.l 
2004 iKitavama et alj|2004l) , we adopt gas density 
around SN progenitor, nsN, as being different from 
the hydrogen gas density in our one-zone galaxy 
model, riH ■ This is because the gas density around 
SN progenitor, ngN, is closely related with dust 
destruction process by reverse shocks driven by 
SNe, as shown in Section [2T2l 



3.2. Star formation law 

We expect that the SFR, ^(t), is roughly pro- 
portional to t~ ir , since a representative timescale 
of the dynamics of the gas disk is t c i T . We assume 
that 

/ H2 (i)M H (i) 



*(*) 



^cir(^vir) 



(18) 



where /h 2 is the mass fraction of molecular hydro- 
gen to th e total gas. We sho uld note that obser- 
vationallv iBigiel et al. ( 2008 ) find that FF? is con- 
verted into stars at a constant efficiency in nearby 
spirals and lGnedin et al" I (l2009l) show that the star 
formation recipe in galaxy formation simulation in 
which star formation occurs only in the molecular 
gas can reproduce the observational correlations 
between SFR and the total gas density. 

3.3. Evolution of gas, stars, and metals 

We calculate the time evolutions of the masses 
of hydrogen and helium gases, M gas , stars, M sta r, 
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and metal of element i, M m> i, in the galaxy by 
using the following equations 

dM,„( t ) _ * (i) M^ +m _ 7sN(() 



dt 

dM star (t) 
dt 

dM mA 



M ISM (t) 

- "lcjocta7SN(i) 

M mti (t) , 



dt 



where Mism(*) 



"*(*) 



MiSM(i) 



;7sn(0 
(19) 



M gas (t) +J2i Mn,i(t)i dejecta = 

and m gas and m mi j are the gas 



mass of hydrogen and helium and the metal mass 
of element i in SN ejecta, respectively. The mass 
returning to the ISM, m e j ec t a (m), w gas (m) and 
(^0 throug h a SN with proge n itor m ass, m, 
are taken from Umeda fc Nomotol ( 2002 ) in the 
case of m = 13, 20, 25 and 30 Mq. m gas , m G j 0Cta 
and m m _i are IMF-averaged TO gas (m), w e j e cta('7i) 
and m mj i(m), respectively. Note that metals con- 
sist of not only heavy elements in gas phase but 
also those in dust grains. 

3.4. Chemistry and cooling 

We follow the time evolution of molecular mass 
fraction, /h 2 , ionization degree, x, and gas tem- 
perature, T. We define the molecular fraction of 
hydrogen as 

/h 2 = (20) 
riH 

where 7ih 2 and na are the number densities of 
molecular hydrogen and hydrogen nuclei, respec- 
tively. The molecular fraction is very important in 
our models, because it determines the final cool- 
ing rate of low-metallicity gas. The metal-free gas 
evolution with chemical r eactions and cooling i s 
studied using the mode l bv lTegmark et al.l (1997 ) 



Hutc hings et al 



(2002) and iHirashita fe Ferraral 
(|2002[ ). We summarize chemical reactions con- 
sidered in this paper and their rate coefficients 
(R n ; n = 1, . . . , 11) in Table [I] The e quations are 
based on IHirashita fc Ferraral (|2002n , but we in- 
clude the effect of the dust size distribution on H 2 
formation and the metal-line cooling process. 
The time evolution of the ionization degree is 



described as 
d.r 



dt 



= xfoRlUH - X i?2™H + ri 2 /o, 



(21) 



where fo = 1 — x — /h 2 is the neutral fraction of 
hydrogen. The terms on the right-hand side are 
the rates of collisional ionization, recombination 
and photoionization. Next, the time evolution of 
the molecular fraction is written as 



d/ H2 
dt 



d/ H2 



dt 



+ 



d/ H2 
dt 



dust 



d,/H 2 
dt 



[d/ H2 " 




d/ H2 " 




dt 


+ 

uv 


dt . 


star 



dost 



(22) 



where the terms on the right-hand side are the 
H2 formation rate in gas phase, the H2 formation 
rate on dust grains, the destruction rate in gas 
phase, and the destruction rate by UV photons, 
and the decreasing rate by star formation, respec- 
tively. These terms are given by 



d/H 2 
dt 

d/a 2 " 
dt 

d/H 2 ' 



dt 

d/H s 



dt 



dust 



dest 



UV 



2/oznH(-Rcff,i 4 

2-Rdust^^H/o, 

-/Ha^H^-Keff,, 

-Ti3/h 2 , 



ReS,2), 



foRio + xRn), 



and 



d/ H2 



dt 



-(1-/ H2 )*(t) 



MismW 



(23) 



The effective formation rates of H 2 including the 



effect of destruction rate of H and H are 



Reff, 



R3R4 



/0-R4 + xR 5 +ri 4 /n H ' 



and 



R, 



ff,2 



R5R7 



/oi?7 + xR?, 



ri5/n H ' 

T+ 



(24) 



(25) 



respectively, and the destruction of ITJ due to H 
collision is 



-Rcff,3 = 



R$Rg 



f R 7 + xR s + ri 5 /n H 



(26) 



We will give the dust-to-gas mass ratio, T>, and the 
production rate of molecular hydrogen via dust 
surface reaction, i?dust , in Section 13.51 and reac- 
tion rates of photo-process, T n (n = 12, . . . , 15), in 
Section [3^61 
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At temperature < 10 4 K, the main coolant is 
molecular hydrogen in low-metallicity gas. The 
cooling rate for molecular hydrogen, Ah., , over the 
range 10K <T< 10 4 K is given by (jGalli fe Palla 
19981) 



log 



Ah 2 (T) 



1 nH«H 2 erg cm 3 s _1 



= -103.0 + 97.59Ti og - 48.05T£ 



-10.80T,L - 0.9032T,! 



(27) 



where T log = log 10 (T/K). iGlover fe Abel <l2008h 
have recently given the H2 cooling rates which in- 
clude H — H2 c ollision and H2 — H2 collision path- 
ways, while the lGalli fe Palla (1998) rates include 
only H — H2 collisions. We assume that the lower- 
limit of gas temperature is the CMB temperature. 

At temperature T > 10 4 K, collisional excita- 
tion, An.ce, and (less importantly) ionization of 
atomic hydrogen, Ajj,ci, are more dominant cool- 
ing process t han molecular hydro gen cooling and 
are given by (jrlaiman et al. 1996 ) 



A H ,cc(T) 



n e -nH erg cm 3 s 1 



7.50x10 



-19 



1 



1 + TV 



exp 



(28) 



and 



A H ,ci(T) 



n c -riH erg cm 3 s 1 



4.02x10 



-19 



i + Ti 



exp 



(29) 

respectively, where T5 is gas temperature in units 
of 10 5 K. 

We consider fine-structure cooling by Ci, Cn 
and Oi, which dominates the thermal evolution 
for number density of g as of interest in this pa- 
per (jOmukai et al.ll2005f ). T he related parameters 
of tra nsitions are given in Hollenbach fe McKeei 
(j!989h . 



3.5. Formation of molecular hydrogen on 
dust grains 

The increasing rate of molecular fraction via 
dust surface reaction is estimated as 



d/H 2 
dt 



= 2i?dust^H/o 



dust 



= ^2 f a fj(a)Tra 2 vSda 

(30) 

where v is the mean thermal speed of hydrogen 
and S is the sticking efficiency of hydrogen atoms. 
We assume that the gas follows a Maxwell ian dis- 
tribu tion so that thermal speed is given by (jKriigel 
20081) 

^ (3D 



IT TOJJ 

Here, we define the dust-to-gas mass ratio, T>, as 

P , E H^/i( a ) da (32) 

. Jo 3n H m H 

The reaction rate of the H2 formation on grains, 
-Rdust, can be estimated as 

Jo V 8a Pj J V 3n H m H / 

(33) 

We adopt S = 0.2 for T < 300 K an d S = for 
T > 300 K (jffirashita fc Ferrarall2002l) . 



3.6. Radiative properties 

We follow photo-processes in chemical reaction 
and heating processes in ther mal evolution by us- 
i ng th e models develo ped by Kitavama fc Ikeuchil 
(|2000l ) and improved in lffirashita fc Ferraral (|2002l ) . 
The intrinsic luminosity is assumed to be equal to 
the total luminosi ty of OB stars whose mass is 
larger than 3 M Q (ICoxll2000l) 



-^uv.o(^) 



dm 



3 M 



dt' L(m)(j){m)y(t-t') 



(34) 

where L(m) is the stellar luminosity as a func- 
tion of stellar mass m. For L(m), we adopt the 
m odel of ze r o-met allicity stars without mass loss 
in Schaerer (|2002l) . We assume the spectrum of 



the incident UV radiation from stars is a power 
law with an index a : 



(35) 



Ajv(^) = -fo^Hl) 



where v is the frequency of photons and io^m) 
is the intensity at the ionization threshold fre- 
quency of neutral hydrogen v^a = 3.3 x 10 15 Hz. 
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Table 1: Reaction rates needed to cal culate th e abun da nce of H2 . The unit of the gas temperature T is K 
unless otherwise stated. Reference: 1. lOmukai (|2000h : 2. lGalli fe PaUal (|l998h 



icnF" s 



No. 
1 



Reaction 



H+ + 2e- 



Rate 

cxp[-32.71 + 13.54 ln(T(eV)) 

- 5.739(ln(T(eV))) 2 + 1.563(ln(T(eV))) 3 

- 0.2877(ln(T(eV))) 4 + 3.483 x KT 2 (ln(T(eV))) 5 

- 2.632 x l(H(ln(T(eV))) 6 
+ 1.120 x 10^ 4 (ln(T(eV))) 7 

- 2.039 x 10- 6 (ln(T(eV))) 8 ] 
exp[-28.61 - 0.7241(ln(T(eV))) 

- 2.026 x 10- 2 (ln(T(eV))) 2 
10- 3 (ln(T(eV))) 3 
10- 4 (ln(T(eV))) 4 
10- 5 (ln(T(eV))) 5 
10- 6 (ln(T(eV))) 6 
10- 7 (ln(T(eV))) 7 
10- 8 (ln(T(eV))) 8 

- 3.071 x 10- 9 (ln(T(eV))) 9 ] 

1.4 x 10- 18 T°- 928 exp(-T/1.62 x 10 4 ) 
4.0 x io- 9 T~ 017 (T > 300); 

1.5 x 10" 9 (T < 300) 

5.7 x 10- 6 T- 1 / 2 + 6.3 x 10~ 8 

- 9.2 x lO^T 1 / 2 + 4.4 x 10~ 13 T 
dex[-19.38- 1.523 log 10 T 

+ 1.118(log 10 T) 2 - 0.1269(log 10 T) 3 ] 

6.4 x 10" 10 

2.0 x lO^T- 1 / 2 
3.0 x 10- 10 exp(- 

1.5 x 10- 10 exp(- 



Ref. 



H 



H+ + e" 



H + 7 



3 


H-f 


e — 


> H- + 


7 


4 


H" 


+ H- 


^H 2 + 


e 


5 


H~ 


+ H+ 


— > 2H 




6 




-H+ - 




-7 


7 




+ H - 


^H 2 + 


H+ 


8 


HJ 


+ e" - 


-> 2H 




9 


H 2 


+ H+ - 




+ H 


10 


H 2 


+ H — 


■» 3H 
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H 2 


■f e" - 


-> 2H4 


- e~ 



2.381 x 
3.213 x 
1.422 x 
4.989 x 
5.756 x 
1.857 x 



1.1- 



a K 



-21050/T) (T < 10 4 ) 
-14000/T) (T > 10 4 ) 



k L = 1.12 x 10~ lu exp(-7.035 x 10 4 /T) 
k H = 6.5 x lO" 7 ! 1 - 1 / 2 

x exp(-5.2 x 10 4 /T)[1 - exp(-6000/T)] 
a = 4.0- 0.416 log 10 (T/10 4 ) - 0.327(log 10 (T/10 4 )) 2 
4.4 x 10- 10 T - 35 exp(-1.02 x 10 5 /T) 



dust 



H + H + grain 



H 2 



gram 



see Section [37J 
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In this paper, we simply s et a = 5 according to 
Hirashita fc Ferraral (|2002h . The normalization of 
the intensity is determined by 



L uv ,oexp(-T disk ) 



I\sv{v)&v, (36) 



where ^ m in is the minimum frequency where 
OB stars dominate the radiative energy of star- 
forming galaxies, and Tdisk is the typical dust 
optical depth in the disk. We assume that 
Vmin = 10 15 Hz. This typical optical depth can be 
simply estimated as 



7"disk — ^disk 



3 



ira 2 fj(a)da (37) 



by assuming that extinction efficiency of dust 
is unity in UV. We calcula te L 12 , fu and Ijr 
from Equation (A20) of iKitavama fc Ikeuchil 
(2000 ) and heating ra t e from Equation (A21) 
of Kitavama fe Ikeuchil ( 2000l ). We summarize 
the cross section for the photo-process in Ta- 
ble [2] The H2 photod is sociat ion cross-section 
is given by lAbel et all (|l997f ). However, if 
the H 2 column density becomes larger than 
10 14 cm" 2 , self-shielding effects become important 
(|Draine fe Bertoldilll996l ). T herefore, H 2 dissoci- 
ation rate, Ti3, is given by ([Hirashita fc Ferrara 



2002) 



T 13 = (4tt)1.1 x 10 8 / uv (3.1 x 10 ib Hz) 
/ f \ -0.75 

X { 1014 cm- 2 J S ' (38) 

where Iuv(3.1 x 10 15 Hz) is in the Lyman- Werner 
band. We should note that adoption of r^isk in 
Equation (j38j) gives an extreme upper bound to 
the self-shielding, so we may overestimate self- 
shielding to internal Lyman- Werner photons by 
H 2 . In this paper, we focus on the effects of dust 
on the protogalaxy, so for simplicity, we set an 
extreme upper bound to the self-shielding. 

We do not consider the Lyman- Werner back- 
ground, since we concentrate in evolution of 
atomic line cooling halos (M v ; r > 10 8 Mq in 
z < 10) in which destruction of molecular hy- 
drogen by the Lyman- Werne r back g round is less 
efficient JCTShea fc Norrnar] 120081 [Susal 120081: 
Wise fc Abelll2008t IWise fc Ceijl2009h ~However. 



10 8 — 10 9 Mq halos are not immune to the Lyman- 
Werner background, just self-shielded at their very 



centers. Not all the baryons will be protected from 
external photodissociating flux and this will affect 
H2 production on dust outside the center of halo. 
In the lower mass halos, the Lyman- Werner back- 
ground may be effective to d i ssociate H 2 molecule 
(iMachacek et al.ll200ll 120031; lYoshida et all 120031; 



Susa 



2008). 



2007tlWise fc AbefaOTUO'Shea fc Normanl 



4. Results 

Our fiducial model assumes M v ; r = 10 9 M© and 
z v i r = 10, and includes the dust destruction model 
by both reverse shocks and forward shocks with 
the ISM density around the SN progenitor, tisn = 
1 cm" 3 (see Table[3]for a summary of our models). 
This dark matter halo forms from a 2.5 a density 
fluctuation. We stop the calculation at z = 5. It 
corresponds to the galaxy age of ~ 0.8 Gyr that is 
before SN II to be dominant source of dust grains. 
The initial mass of gas is Af gas = 1.7 x 10 8 Mq and 
the dynamical timescale of circular motion of the 
gas disk is t c i r (z V n- = 10) = 9.5 x 10 7 yr. Note that 
galaxies with M v ; r ~ 10 9 Mq play a critical role 
in the cosmic reionization, since in the relevant 
redshift range for cosmological reionization, z = 
6 — 15, most of reionization radiation is expected 
to come from galaxies with masses less than ~ 
1O 9 - 5 M (IWise fc Cenll2009h . 



4.1. The dust destruction 

We first show the evolution of a galaxy with 
M v j r = 10 9 Mq and z v i T = 10 for various dust de- 
struction models. To clarify the dust destruction 
effects on galaxy evolution, we first show the result 
of the model without reverse shocks and forward 
shocks (model Clm9), then compare the results 
of the models with only forward shocks (model 
Blm9) and with both forward and reverse shocks 
(model Alm9) to the model without both shocks 
(model Clm9). 

In Figure [5l we show the evolution of vari- 
ous quantities without dust destruction (model 
Clm9). The figure shows the time evolution 
of the molecular fraction, /h 2 , the SFR in 
units of Mq yr" 1 , ^, the stellar mass fraction, 
M s tar/(.Mg as + M star ), the metallicity, Z, in units 
of Zq, total dust cross-section per unit volume, 
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Table 2: Cross-sectio ns for photoionizatio n and p ho todissociation proc ess, where the frequency, v, is in units 
of Hz. Reference: 1, iKitavama k Ikeuchil (l2000f ); 2, lAbel et al.l (|l997t ); 3, iTegmark etafl (| 19971 ) 



TOT 



Reaction 



"RcT 



cross section 

(cm 2 ) „_ 

6.30 x 10- 1S (W3.3 x 10 ib )- 3U 
see equation ([38| 
3.486 x 10- 16 (a; - l) 3 / 2 /x 311 
(x = v/1.8 x 10 14 ) 
7.401 x 10~ 18 

dex(-a; 2 - 0.0302a; 3 - 0.0158a; 4 ) 
(x = 2.762 ln(i//2.7 x 10 15 ) 



v range 
(Hz) 

v > 3.3 x 10 lb 
v > 1.8 x 10 14 
!/ > 6.4 x 10 14 



12 
13 
14 

15 



H + 7 - 
H 2 + 7 - 
H- +7 

H++ 7 



H+ + c 
H + e" 



2H 



H + H 4 



Table 3: the dust destruction model and the main parameters. 



model 


dust destruction 


nsN(cm 3 ) 


Af vir (M Q ) 


Clm9 


no destruction 




10 y 


Blm9 


forward 


shocks 


1 


10 9 


Alm9 (fiducial) 


forward 


and 


reverse shocks 


1 


10 9 


A0.03m9 


forward 


and 


reverse shocks 


0.03 


10 9 


A0.1m9 


forward 


and 


reverse shocks 


0.1 


10 9 


A0.3m9 


forward 


and 


reverse shocks 


0.3 


10 9 


A3m9 


forward 


and 


reverse shocks 


3 


10 9 


Al0m9 


forward 


and 


reverse shocks 


10 


10 9 


A0.1m8 


forward 


and 


reverse shocks 


0.1 


10 8 


AO.lmlO 


forward 


and 


reverse shocks 


0.1 


10 10 


AO.lmll 


forward 


and 


reverse shocks 


0.1 


10 11 


Alm8 


forward 


and 


reverse shocks 


1 


10 8 


AlmlO 


forward 


and 


reverse shocks 


1 


10 10 


Almll 


forward 


and 


reverse shocks 


1 


10 11 


A10m8 


forward 


and 


reverse shocks 


10 


10 8 


AlOmlO 


forward 


and 


reverse shocks 


10 


10 10 


AlOmll 


forward 


and 


reverse shocks 


10 


10 11 



14 
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0.01 - 



0.001 



1e-05 
1e+07 



>p 

M sta /(M aas +M ster ) 

D d 2 








y i* 


— 


J? . 


.#* .1'*** 

,* .!»' 

,.•** # 




,,•«** ^ , 





1e+08 
yr 



1e+09 



Fig. 5. — Time evolution for model Clm9 in which 
dust destruction by the reverse shocks and the for- 
ward shocks are not considered. We show molecu- 
lar fraction, /h 2 , SFR in units of M Q yr -1 , vff, star- 
to-gas mass ratio, M sta , T /(M gas + M staT ), metallic- 
ity, Z, in units of Z & , total dust cross-section per 
unit volume, era. -20 1 m unit of 10 -20 cm" , and 
dust-to-gas mass ratio, Z?d.-2, in units of 10 -2 . 
The horizontal axis shows the age of galaxy in unit 
of yr from the formation rcdshift, z v - n - = 10. 

cr d.-20i hi unit of 10 -20 cm -1 , and the dust-to-gas 
mass ratio, in units of 10~ 2 . The defini- 

tion of .Dd.-2 is convenient for comparison with 
the MW value of dust-to-gas mass ratio. In the 

x 1(T 2 in 



MW the dust-to -gas mass rat io is 0.5 
the diffuse ISM dDraind 120091) and .9 



x 10" 



molecular clouds (jPollack et al.lll994f ) 



The molecular fraction, /h 2 , is very important, 
since it determines SFR and controls galaxy evo- 
lution. The molecular fraction reaches /h 2 ~ 
1 x 10~ 3 around t ~ 10 7 yr. This results is robust 
for all models, since in this stage H2 formation 
in the gas phase is dominant over that on dust 



grain s (ITegmark et al.l Il997t iHirashita fc Ferrara 
20021 ). The gas temperature rapidly drops be- 
low 200 K before 10 7 yr. Then, the molecular 
fraction rapidly increases from t ~ 10 8 yr and 
reaches ~ 0.83 at the galaxy age of ~ 0.8 Gyr 
(z = 5). This is due to the enhancement of H2 
formation on dust grains by increase of <7d,_20- 
For t > 3 x 10 7 yr, a d ^ 2 o > 1-1 x 10~ 4 'and 



[d/H 2 /dt] dust exceeds [d/ H2 /df] star . For <r d ^ 2 o > 
0.001, the increase of molecular fraction enhances 
the star formation. The cycle of the H2 formation 
on dust, the star formation and the dust forma- 
tion by SNe, significantly accelerates galaxy evo- 
lution, such as rapid increase of the stellar mass 
fraction, M star /(M sas + M star ). At the galaxy 
age ~ 0.8 x 10 9 yr, the stellar mass fraction goes 
up to M star /(M ga8 + M star ) ~ 0.60. The SFR, 
&(t), decreases from the time when M star /(M gas + 
Mgtax) ~ 0.45, since gas mass decreases signifi- 
cantly. The active star formation causes the for- 
mation of dust grains and metals. At t ~ 0.8 Gyr, 
the total dust cross-section, (Td,-20, goes up to 2.3 
and the metallicity, Z, goes up to 3.3 x 10 _1 Z®. 
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Fig. 6. — Same as in Figurc[5]but for model Blm9 
in which the dust destruction only by the forward 
shocks is considered. 

Figure [6] shows the results of the galaxy model 
with the dust destruction by only the forward 
shocks (model Blm9) to illustrate the effects of 
dust destruction by forward shocks on the galaxy 
evolution. In this model tt-sn = 1.0 cm~ 3 . The 
dust destruction by forward shocks slightly af- 
fects the dust-to-gas mass ratio, Aj,-2i after the 
galaxy age of ~ 5 x 10 s yr. This is because the 
destruction by forward shocks is roughly propor- 
tional to the dust-to-gas mass ratio (see Equation 
©), and dust grains are destroyed significantly 
for -Dd,-2 ^ 0.1 in this case. The SFR decreases 
from the time when M stair /(M g&s + M star ) ~ 0.4. 
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At ~ 0.8 Gyr, the molecular fraction reaches 
/h 2 ~ 0.51 and the stellar mass fraction reaches 
M sta r/(A/ gas + M star ) ~ 0.47. 
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Fig. 7. — Same as in Figurc[5]but for model Alm9, 
in which the destruction by reverse and forward 
shocks in the case of Usn = 1 cm~ 3 
is considered. 

In Figure we show the results of our fiducial 
model Alm9 in which we include the dust destruc- 
tion by both reverse and forward shocks in the case 
of nsN = 1 cm" 3 . The molecular fraction reaches 
to /h 2 ~ 1 X 10~ 3 around t ~ 10 7 yr. This is 
similar to Clm9 and Blm9. After t > 10 7 yr, the 
molecular fraction evolution is quite different from 
models Clm9 and Blm9. The molecular fraction 
declines slowly until t ~ 2 x 10 8 yr. After ~ 
2 x 10 8 yr, the molecular fraction increases slowly 
with increase of the dust mass. This is due to the 
H2 formation on the dust grains. The molecular 
fraction reaches only ~ 2.0 x 10~ 3 at t ~ 0.8 Gyr. 
This is because dust destruction by reverse shocks 
is very effective and hence results in suppression 
of H2 formation on dust grains. On the other 
hand, forward shocks hardly affect the evolution 
of dust size and dust mass, since the destruction 
of forward shocks can change dust mass only for 
large dust-to-gas mass ratio, Ai,-2 ^ 10 _1 . At 
t ~ 0.8 Gyr, the stellar mass fraction reaches only 
Af st ar/(M g as + M s tar) ~ 4.5 x 10" 3 , which is much 
less than the model without reverse shocks shown 
in Figure [5] (model Clm9) and Figure |TJ] (model 



Blm9). 

We illustrate the difference in the H2 forma- 
tion rate among models Clm9, Blm9, and Alm9 
as follows. The H2 formation rate depends not 
only on the total dust mass but also on the dust 
size distribution. In models Clm9 and Blm9, 
the dust mass prod uced b y a SN II without re- 
verse shock is J2j f-M% -(a)da = 0.48 M Q and 
in the Alm9, the dust mass injection into ISM 
through a rev erse shock with ngN = 1-0 cm~ 3 
is Ej/ Ml°(a)da = 0.15 M Q . The ratio of 
the mean dust area to the mean dust volume is 
(a 2 ) /(a 3 ) = 1.5 x 10 5 cm -1 before the reverse 
shock destruction. After the reverse shock de- 
struction, (a 2 )/(a 3 ) = 4.2 x 10 4 cm" 1 . This is 
a measure of dust area per the dust volume and 
is also a measure of H2 formation rate of the dust 
surface. Small (a 2 ) /(a 3 ) leads to a low H2 for- 
mation rate. This is the reason why model Alm9 
shows smaller H2 fraction than Blm9 and Clm9. 
The dust destruction by reverse shocks changes 
not only the dust mass but also the grain size dis- 
tribution, and as a result drastically suppresses 
star formation in the galaxy. 

In Figure [5J we show the evolution of H2 for- 
mation rate in gas phase, [d/H 2 /di] gas , the H2 
formation rate on dust grain, [d/H 2 /di]dust) the 
H 2 destruction rate by UV photons, [d/H 2 /dt]uV: 
and the H2 decreasing rate by the star formation, 
[d/H 2 /di] s tar, normalized to the total formation 
rate, [d/ H2 /dt] g as + [d/H 2 /dt] dus t, in the model 
Alm9. At t — 4 x 10 7 yr, the H 2 formation 
rate on dust grains exceeds the rate in gas phase. 
However, the H2 formation rate on dust grain 
is less than the H2 decreasing rate by star for- 
mation at this epoch. At t ~ 1.6 x 10 8 yr, the 
H2 formation rate on dust grains exceeds the H2 
decreasing rate by star formation. The forma- 
tion on dust grain becomes the dominant process 
among all of H2 formation and destruction pro- 
cesses. From this time when <7d.-20 ^ 5.2 x 10 -5 , 
molecular fraction, /h 2 , starts to increase. The 
molecular destruction rate by UV photons does 
not exceed the molecular decreasing rate by 
the star formation after t ~ 4 x 10 7 yr. Af- 
ter t - 4 x 10 7 yr, Juv(3.1 x 10 15 Hz) ~ 
2.0 - 4.8 x 10~ 20 erg s" 1 cirT 2 Hz -1 ster" 1 , 
which corresponds to J21 = 20 — 40 where 
J21 is in units of Itjv(3-1 x 10 15 Hz) = J 2 i x 

1Q-21 crg g-1 cm ~2 jj z -1 ster -l_ fj 0tc t h at J 21 
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Fig. 8. — Formation and destruction rates of H2 
in model Alm9. The evolution of molecular for- 
mation rate in gas phase (gas), the molecular for- 
mation rate on dust grain (dust), the molecular 
destruction rate by UV photons (UV), and the 
molecular decreasing rate by the star formation 
(star), normalized to the total formation rate, are 
shown. The molecular formation on dust grains 
becomes dominant, compared with the other chan- 
nels after t ~ 1.6 x 10 8 yr. The molecular fraction, 
/h 2 , and total dust cross-section per unit volume 
in unit of 10 -20 cm , <7d,-20) are the same as in 
Figure [7] 



in this model is higher than the values, J 21 ~ 1, 
for the Lyman- Werner background in the redshift 
of 5 < z < 10 sugg ested in recent papers (e.g. 
Greif fc Brorrrm]|2006l) . The destruction process in 
gas phase does not affect the evolution of molecu- 
lar fraction significantly after t > 1 x 10 7 yr. 



4.2. The ISM density around SN 

The dependence of time evolution of SFR on the 
ISM density around SN is presented in Figure [9] 
We can see that after t ~ 5 x 10 7 yr higher density 
around SNe progenitors results in lower molecular 
fraction, and hence lower star formation efficiency. 
The SFR, ^(t), is independent of rigN before t ~ 
5 x 10 yr, because H 2 forms predominantly in the 
gas phase. In the model without reverse shock 
destruction, the SFR increases from ~ 10® yr and 
saturates around 5 x 10 s yr. This is because the 
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Fig. 9. — Time evolution of SFR in unit of 
Mq yr -1 in the model with dust destruction for 
various nsN = 0.03, 0.1, 0.3, 1.0, 3.0 and 10 cnr 3 . 
The values of tisn are used in the panel as in Fig- 
ure [TJ The model without dust destruction by the 
reverse shocks is also shown. The horizontal axis 
shows the age of galaxy in unit of yr from the for- 
mation redshift, z v w = 10. 
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Fig. 10. — Molecular faction evolution, /h 2 , for 
the model of different densities around a SN II, 
"ism.sn = 0.03, 0.1, 0.3, 1.0, 3.0 and 10 cm" 3 , 
which are used in the panel as in Figure [I] and 
the model without reverse shocks. The horizontal 
axis shows the metallicity, Z, in unit of Z Q . 



gas is consumed by the star formation. In model 
A10m9 (ngN = 10 cm -3 ), SFR is suppressed until 
t ~ 0.8 X 10 9 yr. 

We should note that it is probably that tt-sn < 
1 cm -3 for Pop III stars in the mass range of 
20 — 40 M Q , since Pop III stars are massive and 
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Fig. 11. — Molecular faction evolution, /h 2 , for 
the model of different densities around a SN II, 
msM.SN = 0.03, 0.1, 0.3, 1.0, 3.0 and 10 cm" 3 , 
which are used in the panel as in Figure [TJ and 
the model without reverse shocks. The horizontal 
axis shows the total dust cross-section per unit 
volume, (Td,-20 5 m unit of 10~ 20 cm -1 . 



can photoevaporate the clouds in which they form, 
and that ionized flows evacuate the dense gas 
around the stars to well below nsN = 1 cm -3 
(|Whalen et all l2004t iKitavama et all |2004 . In 



this case, considering circumstellar densities of 
10 cm -3 and greater is not relevant to dust evolu- 
tion in the SN remnant. In this paper, we consider 
"■sn > 5 cm -3 for completeness. If the stars are 
forming at lower redshift and are enriched, they 
will have stellar winds that also sweep away cir- 
cumstellar gas to low densities. 

We note that in usual star formation recipe in 
both numerical simulations and analytic models, 
SFR is assumed to increase with the local gas den- 
sity and our results show that the SFR is strongly 
affected by ngN- We will discuss the effects of tisn 
on the SFR in more detail in Section [5j 

In Figure QUI we show the change of the molec- 
ular fraction, /h 2 , with metallicity, Z, for various 
ngN- In usual chemical evolution models, Z is a 
key indicator of the galaxy evolution. However, as 
shown in this figure, /h 2 does not solely depend 
on the metallicity. For Z > 5 x 10~ 4 Zq, /h 2 is 
large in models with small tisn- This is because 
-Darea/-Mmctai is large (small) in models with small 
(large) nsN for the same Z as shown in Figure |3l 

The molecular fraction is well described by 



the total dust cross-section per unit volume for 
°d.-20 <; 0.001. In Figure HU we show the 
evolution of molecular fraction in terms of total 
dust cross-section per unit volume. For era, -20 > 
1 x 10~ 3 , H2 formation on dust grains dominates 
/h 2 evolution as shown by the convergence of all 
models. 



4.3. The dark matter halo mass 
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Fig. 12. — The stellar mass fraction in the galaxy, 
Af s tar/(Afgas+A/star), for the various virial masses, 
M v ir of the models with nsN = 0.1 cm -3 (solid 
line, model A0.1m8, A0.1m9 , AO.lmlO, and 
AO.lmll), with tiism.sn = 1-0 cm~ 3 (dashed line, 
model Alm8, Alm9, AlmlO, and Almll), and 
with nisM.SN = 10 cm~ 3 (dashed line, model 
Al0m8, Al0m9, AlOmlO, and AlOmlO). 

Finally, we show the stellar mass fraction, 
AW(A/ gas + AW) for Af vir = 10 8 , 10 9 , 10 10 
and 10 11 Af at z vir = 10, in Figure [TJ We 
note that dark halos of virial masses, Af vlr = 10 8 , 
10 9 , 10 10 and 10 11 Af Q correspond to the density 
fluctuation of 2. Oct, 2.5ct, 3. Oct and 4.1ct, respec- 
tively. Af s tar /(Af gas + Af sta r) is large for large 
M v j r . This is explained as follows; the gas cools 
to the CMB temperature ~ 30 K in all M v - lr after 
t ~ a few x 10 7 yr although T VK increases with 
Af v i r , so that the final gas density becomes higher 
because of smaller if/r^sk in larger M v i r (see 
Equation (|15[) V This results in more rapid molec- 
ular formation in a larger Af v i r halo (see Equation 
([23)) ). The rapid molecular formation enhances 
the star formation and as a results, causes the 
large stellar mass fraction. We should note that 
in higher z vir , Af star /(Af gas + Af star ) is larger for 
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the galaxies with T vu > 10 4 K, since rotation 
timescale become smaller in our model. 

5. Discussion and conclusions 

We have investigated the dust size evolution 
and the resulting H2 formation on dust grains in 
the galaxies in the early Universe (z > 5). Our 
dust evolution model includes the dust production 
by SNe II and the effects of dust destruction by the 
reverse shocks and forward shocks driven by SNe. 
In the galaxy model, we follow the chemical net- 
work of H 2 formation and the thermal evolution 
of gas taking into account cooling by H 2 , atomic 
hydrogen and metals, and heating by the stellar 
radiation. The SFR is assumed to be proportional 
to the mass of molecular hydrogen and the rota- 
tion timescale of galactic gas disk. The evolution 
of size distribution of dust has been taken into ac- 
count for the first time in this paper to investigate 
the influence on star formation activity in galaxies. 

We obtain three main results. First, we show 
that H2 formation is suppressed by the dust 
destruction, especially by the reverse shocks in 
SNRs. Resulting molecular fraction in the galaxy 
model with dust destruction by both reverse 
shocks and forward shocks (see Figure [7]) is 2.5 
orders of magnitude less than the galaxy model 
without both destruction processes (see Figure [5]) 
and is 2 orders less than the galaxy model con- 
sidering dust destruction by only forward shocks 
(see Figure O for tisn = 1-0 cnr 3 . Our results 
show that dust size evolution has great effects on 
the early galaxy evolution. 

The dust destruction by reverse shocks is essen- 
tial in our galaxy model and has more dominant 
role than forward shocks, since the reverse shock 
destroys small dust grains earlier than the forward 
shock. Dust destruction by a reverse shock be- 
comes more (less) prominent for a higher (lower) 
nsN- In early galaxy evolution, considering dust 
destruction by reverse shocks is very important to 
H2 formation. We note that forward shocks affect 
dust size distribution in large dust-to-gas mass ra- 
tio, T> > 10 -3 , as shown in Section |4~T1 

Second, we show that the SFR strongly depends 
on the ISM density around SNe progenitors, tisn, 
which determines the efficiency of dust destruction 
by reverse shocks. As noted in Section |3l we treat 
??sn different from the mean gas density in our one- 



zone galaxy model, tih , taking into account the ef- 
fect of photo-ionization by SNe progenitors. The 
molecular fraction is different among the models 
with various nsN> even for the same Z (see Fig- 
ure [TUJ) , and it solely depends on total dust cross- 
section per volume from the time when H2 for- 
mation on dust grain becomes the most dominant 
process than the other processes of H2 formation 
and destruction (see Figure fTTj) - The ratio of total 
dust cross-section to total metal mass presented in 
Figure [3] is very useful for a numerical simulation 
of galaxy formation with dust size evolution. 

A higher ISM density around SN progenitors, 
nsN results in lower SFR in the early galaxy evo- 
lution. In conventional galaxy formation theory, 
SFR has been assumed as an increasing function 
of gas density. Since in our one-zone model, we 
simply treat ngN as a parameter, it is very inter- 
esting to resolve the ionized region around a SN 
progenitor by high-resolution radiation hydrody- 
namic simulation taking into account the effects 
of ionization heating by massive SNe progenitors. 

Finally, we show that galaxy evolution depends 
on M v i r and show that the stellar mass fraction, 
M s t ar /(M gas + Af s tai), is monotonically increas- 
ing functions of M V - 1T . The halo number density 
in the redshift range 5 < z < 10 is a decreas- 
ing function of the halo virial mass, 10 s M Q < 
M v j r < 10 11 M Q . In most papers on the reion- 
ization, M star /(M gas + M star ) is assumed to be 
constant for various halo mass, M v ; r . The differ- 
ence of stellar mass fraction with different M v j r 
in our results is important for galaxy evolution 
in the early Universe and should be considered 
in the cosmic reionization process. We will study 
this effects taking into account Population III star 
formation and the H2 dissociation by the Lyman- 
Werner background in a forthcoming paper, since 
strong Lyman- Werner background that is effective 
to dissociate molecular hydrogen, to delay gas con- 
densation, and to suppress SF especially in the 
lower mass H2 cooling halos, is expected in the 
cosmic reionization era (Yamasawa et al. in prepa- 
ration). 



Dust size distribution governs the extinc- 
tion curve and reemitting IR spectrum (e.g. 



Silva & Denese 


1998; Granato et al.l2000l Takeuchi et al. 


2005; Li et al.l 


2008). A ground-based intcrfero- 



metric facility, the Atacama Large Millimeter Ar- 



ray (ALMA) can be used to study high-rcdshift 
galaxies, since redshifted dust emission can be ob- 
served with ALMA. We can study dust emission 
and resulting SED of our galaxy model using our 
results of the dust size distribution and test our 
galaxy model by comparing to observational data 
in the future. 

Our galaxy model can be applied to the lower 
redshift (z < 5) star formation history, taking 
into account dust production by low-mass stars 
and g rowth of dust due to accr etion of metal in 
ISM (jlnoue I [200I iDraind liooj ) . The accretion 



timescale, r acc , depends on the metallicity, and is 
given by r arc = T acc , Z e /Z, where T acC;0 ~ 10 8 yr 



( Inoue 20031 ). If we assume Z ~ 10~ 2 Zq corre- 
sponding r acc ~ 10 10 yr, then the accretion time 
is too long to affect dust size distribution for a 
cosmic time < 1.2 Gyr (z > 5). At a galaxy age 
of > 1 Gyr, dust production by low-mass stars 
should affect the dust size evolution, because after 
that epoch dust is supplied from low-mass stars 
as well as SNe II. We will include such processes 
to investigate low-redshift galaxies in the future 
work. 

Cosmological simulation of galaxy formation in- 
cluding our dust formation and evolution model is 
needed, since it is widely understood that most 
z ~ 10 galaxies were not clear disk galaxies 
(I Johnson et al.ll2008t iGreif et alJlioiot IWise et all 



2010l ). We will study first galaxy formation by 



cosmological simulation including our dust model. 
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